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Symulacje komputerow 
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Gra w życie 


e Najbardziej znany automat komórkowy wymyślony 
w 1970 przez Johna Conwaya 


° Dwuwymiarowa nieskończona plansza (każda komórka 
ma ośmiu sąsiadów ) z dwustanowymi polami (martwa 
lub żywa) 

e Reguły (zmiana stanu pola zależy od ilości jej sąsiadów): 


1. W martwym polu rodzi się żywa komórka, jeżeli ma 
dokładnie trzech żywych sąsiadów 

2. Warunkiem pozostania przy Życiu jest posiadanie 
dwóch lub trzech żywych sąsiadów; w przeciwnym 
przypadku umiera z samotności lub „zatłoczenia” 


Gra w życie 


° Do pokazów będę używał implementacji 
Iwony Białynickiej-Biruli z książki 
Modelowanie rzeczywistości 
napisanej przez pisia 41101-1 + ТЕЛЕТ 
Iwo Białynickiego-Birulę 
(WNT, 2007) 


«"E Modeling Reality - Conway 


Starsza wersja dostępna na WWW: 
http://www.wiw.pl/modelowanie/conway.asp 


Gra w życie 


Analiza reguł na przykładzie stanu złożonego z 3 pól 
(określamy ilość sąsiadów) 


W efekcie powstaje układ oscylujący u 


NEM 
między dwoma stanami (okres = 2): Е 


Gra w życie 


° Stany stałe (stabilne, niezmienne) 
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blok koniczynka staw kryształ łódź bochenek 


° Oscylatory 


blinker (2) żabka (2) fontanna (15) krokodyl (16) 


Gra w życie 


e Rozwój układu może kończyć się 
stanem stałym, oscylatorem lub śmiercią kolonii 


M = MH M 
KBE = m а 
MER m= KIER wszystkie te stany ewoluują 
E ШЕ NEM do takiego samego oscylatora 
NEM K NEM 


złożonego z czterech blinkerów 


E nieuchronna śmierć kolonii 


ewolucja tych stanów kończy się kryształem 


Gra w życie 


° Przykład dłuższej ewolucji z oscylatorem na końcu 
(Delta 02/1977) 


° Przykład ewolucji n" 
do stałego stanu n" 
(fabryka bloków, si 


żniwa = harvest) =" 


Gra w życie 


° Generowanie stałego układu obejmującego całą 
powierzchnię 


Gra w życie 


e Układy oscylujące z przesunięciem (statki) 


Szybowiec (glider) 


Dakota 


Inne układy: szalupa, fabryki szybowców (działa), big ship, puffery, 
R-Pentomino 


Gra w życie 


° Stan komórki jest deterministyczny: jednoznacznie 
zależy od ilości jej sąsiadów 
(nie zależy od ich ułożenia) 

« Jednak ewolucja w tym automacie 
nie jest odwracalna tzn. wiele stanów może 
prowadzić do tego samego końca (zob. oscylatory) 

° [stnieją również stany, których nie można uzyskać w 
efekcie ewolucji żadnego możliwego stanu 


początkowego (tzw. rajskie ogrody) 


Automaty komórkowe 


° Ogólna definicja automatów komórkowych: 

— n-wymiarowa siatka komórek o; (i — zbiór indeksów) 

— komórki mogą przyjmować stany z zadanego przestrzeń 
stanów (żywy/martwy, kolory, liczby naturalne, itp) 

— reguły określające w jaki sposób stan komórki w kolejnej 
chwili czasu zależy od stanu całego układu w chwili 
poprzedniej o; (t +1) = F(t (0) }, j należy do otoczenia i) 

— jeżeli ewolucja zależy od zmiennej losowej, automat 
nazywany jest probabilistycznym 


Gra w życie: п = 2, zbiór stanów to +0, 1} 

reguła: F = 1 dla trzech żywych sąsiadów; 

dla dwóch żywych sąsiadów komórka zachowuje stan F = J; 
w pozostałych przypadkach F = 0 (umiera) 


Automaty komórkowe 


° Otoczeniem komórki i nie muszą być wszystkie 
komórki z najbliższego sąsiedztwa 


Sąsiedztwo Moore'a i von Neumanna 
dla promienia równego jednej komórce 


Automaty komórkowe 


° W symulacjach fizycznych (materiały sypkie, gazy, 
ciecze) stosuje się sąsiedztwo Margolusa 1 reguły 
ustalające zachowanie jednocześnie czterech komórek 


° Bloki komórek przesuwają się w lewo 1 w prawo 


° То wymusza parzyste 1 nieparzyste kroki ewolucji 


Automaty komórkowe 


° [nne znane automaty komórkowe dwuwymiarowe: 


— mrówka Langtona - dwuwymiarowa maszyna Turinga, 
której program jest następujący: 
wyróżniona komórka (mrówka), która posiada kolor biały lub 
czarny i kierunek (N, E, S, W), obraca się w lewo jeżeli 
znajduje się na polu białym, a w prawo, jeżeli na czarnym. W 
każdym przypadku zmienia kolor na przeciwny 1 wykonuje 
ruch zgodnie z kierunkiem 


Automaty komórkowe 


° [nne znane automaty komórkowe dwuwymiarowe: 


— Model Greenberga i Hastingsa 
automat samoporządkujący się, który z losowego ułożenia 
stanu komórek (możliwe stany to +0, 1,2} ewoluuje po 
odpowiednio dużej liczbie kroków do stanu oscylującego 


— Modeluje reakcję Biełorusowa-Żaboryńskiego (chemia) 
zmiana barwy trzech zmieszanych roztworów wodnych 


t=0 1= 58 t= 105 t= 155 1= 205 
t= 255 1= 305 1= 355 1= 405 t= 455 


Aplet Java (psoup.math.wisc.edu/iava/iqh.html#T 


Automaty komórkowe 


° [nne znane automaty komórkowe dwuwymiarowe: 


— model Greenberga 1 Hastingsa 


Reguły: 

komórka w stanie 2 (stan przesycenia) przejdzie do stanu 
0 (równowagi) bez względu na swoje otoczenie 

(von Neumanna) i nie wpływając na nie 

komórka w stanie 1 (stan przejściowy) zmienia się w 2 

1 dyfunduje na swoje otoczenie tzn. zmienia stany О w 
otoczeniu na stany 1 

komórka w stanie 0 (stan równowagi) pozostają w stanie o 
ile nie są zmieniane na mocy poprzedniej reguły 


Automaty komórkowe 


e Najlepiej zbadane są automaty jednowymiarowe 


° Słynna książka Stephena Wolframa (auto 
programu Mathematica 1 witryny MathWorld) 
pt. A New Kind of Science 


e Hipoteza o równoważności 


automatów jednowymiarowych F ZBY 
С KIND OF 


Z MASZYDAMIE гапира (С SCIENCE 
(o nich za chwilę) н 


Automaty komórkowe 


° Próba wprowadzenia przez Wolframa 
podziału automatów (jednowymiarowych): 

— Klasa I: Automaty niezmienne - ewoluują do stałego 
stanu homogenicznego, np. do śmierci całej kolonii 

— Klasa II: Automaty ewoluujące do stanu stałego lub 
okresowego (oscylatora) 

— Klasa III: Automaty chaotyczne — nieuporządkowane 
lokalnie 1 globalnie (nie mające wzorca zachowania) 

— Klasa IV: Pozostałe automaty ewoluujące według 
ciekawych 1 unikalnych praw 


Automaty komórkowe 


Automaty komórkowe mogą być stosowane do wielu 
zagadnień biologii, chemi 1 fizyki 

Nie udał się jednak zamiar udowodnienia 
równoważności z równaniami różniczkowymi 

W fizyce: symulacje ciał sypkich (piasek), gazów 1 
cieczy, przeciekanie przez ciała porowate (percolation), 
model Isinga, różnego typu przejścia fazowe 

I wszystkie te zagadnienia, w których nie można 
zaniedbać lokalnych korelacji (np. równanie kinetyczne 
Boltzmana 1 przybliżenie chaosu molekularnego) 


Poza tym: pożary lasu, korki uliczne, modele społeczne 


Paradygmat systemowy 


Proste reguły rządzące zachowaniem komórek prowadzą do 
złożonych (1 niespodziewanych) zachowań całych kolonii 


Zaczynając od prostego układu, poprzez chaotyczny rozwój można 
nieraz osiągnąć układ bardzo złożony 


Układy złożone - zbudowane z modułów (komórek) związanych 
sprzężeniami zwrotnym. Dzięki temu własności całości są Inne niż 
poszczególnych elementów. 


„Całość to więcej niż prosta suma części” (Arystoteles) 


W sprzeczności z redukcjonizmem (Kartezjusz), w którym 
wszystko można zrozumieć przez analizę 1 późniejszą syntezę 


Holizm (Smuts, XX wiek) = całość nie redukuje się do części 
Synergia — wzmocnienie dwóch lub więcej czynników układu 


Adaptacja i ochrona układu dzięki nadmiarowości 


Paradygmat systemowy 


Cybernetyka (1948-19535) to nauka o efektywnych systemach 
sterowania (organizacji), mechanizmach kontroli i związanej z tym 
komunikacji między układem sterowanym 1 jego kontrolerem. 

Ma zastosowanie w teorii maszyn 1 zwierząt. 

Ogólna teoria systemów (1950) — ogólny opis organizmów 
żywych, ale również społeczeństw 1 układów sztucznych 

Teoria katastrof (1970) — teoria przejść nieciągłych tzn. takich, w 
których ciągła zmiana parametru kontrolnego zmienia jakościowo 
własności całego systemu. Przejścia fazowe. 


Teoria chaosu (1980) — badanie układów dynamicznych, których 
ewolucja jest bardzo wrażliwa na warunki początkowe. 


Teoria złożoności (1990) 


Źródło: prof. Marian Ostwald (http://www.sms.am.put.poznan.pl/?site=skrypty) 


Paradygmat systemowy 


Teoria złożoności (science of complexity) 


Teoria informacji Cybernetyka Prakseologia 


Ogólna teoria systemów 
Badania systemowe 
Badania operacyjne 


Analiza systemów Inżynieria systemów 


Ogólna teoria 
optymalizacji 


Optymalizacja 


systemów i konstrukcji inżynierskich. 
Optymalizacja działań, wyrobów, usług, 


Źródło: prof. Marian Ostwald (http://www.sms.am.put.poznan.pl/?site=skrypty) 


Paradygmat systemowy 


Krytyczny poziom złożoności powyżej którego układ nabiera 
nowych jakościowo cech (np. inteligencja lub świadomość mózgu) 


1 2 
dividing zygote, 
implantation and 

gastrulation 


external genitalia 


| Ба 
not — 


susceptible to ma 2) | | 
eratogens 
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Automaty komórkowe a fizyka 


Ciała sypkie (piasek) 
a) om "ЫЙ т 
Otoczenie Margolusa 
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"M 
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B. Choppard, M. Droz, Cellular Automata Modelling of Physical Systems, Cambridge University Press 1998 


Automaty komórkowe a fizyka 


Ciała sypkie (piasek) 


Ruchome Piaski 1.0 


Tryb Edycji 
(# Sypanie piaskiem 
С Edycja brzegów 
Czyść Wszystko 


Zapisz 


Wczytaj 


Start | Step Stop | Reset 


Maciej Matyka: http://panoramix.ift.uni.wroc.pl/”maq/pl/automat.php 
Aplet Java: http://schuelaw.whitman.edu/JavaApplets/SandPileApplet/ 


Automaty komórkowe a fizyka 


Symulowanie gazu 


Maciej Matyka: http://panoramix.ift.uni.wroc.pl/”maq/pl/automat.php 


Automaty komórkowe 


° Pomysły na projekty zaliczeniowe 
— [mplementacja modelu Greenerga-Hastingsa 
— Implementacja przesypywanmia piasku 
— [mplementacja rozchodzenia gazów 
— [mplementacja maszyny Turinga 


— Przygotowanie programu MT do sortowania О 1 1 


° Projekty mogą się powtarzać o ile będą 
realizowane w różnych technologiach 
(preferowane aplety Java) 


Dodatek: maszyna Turinga 


° Budowa maszyny Turinga 
— taśma (pamięć trwała), w jej komórkach zapisane są 
wartości z ustalonego alfabetu np. A - Z 1 kropka (,,.”), 
poza literami R 1 L, które będą kodować ruch głowicy 
— głowica odczytująca 1 zapisująca 
stany głowicy, czyli rejestru, wybierany jest 
z osobnego alfabetu np. a - z 


— program (czwórki) (stan 


Dodatek: maszyna Turinga 


e Program to czwórki znaków np. 


bieżący stan głowicy qASs nowy stan rejestru 


(stan rejestru) 


bieżąca wartość nowa wartość zapisana 
(wartość na taśmie na taśmie w tym miejscu 


wskazywana przez 
głowicę) 


e Jeżeli alfabety taśmy 1 głowicy są rozłączne, stan maszyny 
(stan taśmy 1 rejestru oraz pozycję głowicy) można 
jednoznacznie zapisać przez AJAAABAA 
(taśma AAAABAA, rejestr q, pozycja głowicy п = 2) 

° Program nie może mieć dwóch czwórek (linii) 

o takich samych dwóch pierwszych znakach (jednoznaczność) 


Dodatek: maszyna Turinga 


° Prosta maszyna Turinga 
— taśma: AAAABAA 


— głowica: początkowy stan rejestru q, pozycja n = 2 


— program (3 linie): qASS — zamień wartość z A na S i ustaw rejestr na s 


° Output: 
AqAAABAA 
AsSAABAA 
ASqAABAA 


ASSSSqBAA 
ASSSSBbAA 


SSRq — przesuń głowicę w prawo i ustaw rejestr na q 


qBRb — przesuń głowicę w prawo i ustaw rejestr na b 


pasuje qASs 
pasuje qSRq 
pasuje nA Se 


pasuje 


Dodatek: maszyna Turinga 


Koniec działania maszyny — gdy nie ma czwórki 
pasującej do bieżącej wartości 1 stanu rejestru 
Można również wprowadzić wyróżniony stan 
rejestru, który sygnalizuje koniec programu 

(w naszym przykładzie jest to stan b) 


Udowodniono, że zwykłe komputery są 
równoważne maszynom Turinga! 


Studiowanie problemu rozwiązywalności 


Dodatek: maszyna Turinga 


° Filozofia 1 psychologia poznawcza: maszyna 
Turinga jest wygodnym narzędziem przy 
precyzowaniu pojęć 1 problemów procesu 
poznawania 1 sztucznej inteligencji 

° Słynne zagadnienie: czy człowiek jest maszyną 
Turinga 1 jeżeli nie, jak ich odróżnić za pomocą 
skończonej liczby pytań (psychiatra w Emacs) 


° | wiele innych... 


Jeszcze raz 


° Pomysły na projekty zaliczeniowe 
— [mplementacja modelu Greenerga-Hastingsa 
— Implementacja przesypywanmia piasku 
— [mplementacja rozchodzenia gazów 
— [mplementacja maszyny Turinga 


— Przygotowanie programu MT do sortowania 0 1 1 


• Projekty mogą się powtarzać o ile będą 
realizowane w różnych technologiach 
(preferowane aplety Java) 


